New approach to phase transitions in black holes 
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We develop an analogy between fluids and black holes to study phase transitions in the latter. 
The entropy-temperature graph shows the onset of a phase transition without any latent heat. The 
nature of this continuous (higher order) phase transition is examined in details. We find that the 
second order derivatives of the free energy diverge at the critical temperature. Also, the transition is 
smeared instead of sharp, so that the usual Ehrenfest's scheme breaks down. A generalised version 
of this scheme is formulated which is shown to be consistent with the phase transition curves. 



Although black holes are well known thermody- 
namical systems, they are not well understood. In 
fact there is no microscopic or statistical description 
behind their thermodynamical behaviour. Studying 
phase transition in black holes is very important as it 
could manifest the underlying micro-structure of black 
hole thermodynamics. Although there are some inves- 
tigations jj, these are more concerned with the geo- 
metrical rather than the thermodynamical aspects of 
phase transition. 

In this paper we systematically develop a new ap- 
proach to study phase transitions in black holes from 
a thermodynamical viewpoint. To do that we first 
build an analogy between fluids and black holes that 
parallels the analogy between between fluids and mag- 
netic systems. From this analogy we then exploit the 
known results in fluids to infer the corresponding sit- 
uations in black holes. Following this approach we 
find that the entropy-temperature graph (as well as 
the Gibb's energy-temperature graph) of the Kerr- 
AdS black hole shows the onset of phase transition 
without any latent heat. Moreover this transition is 
smeared round the critical temperature. The details 
of this continuous (higher order) phase transition are 
examined by looking at various plots of entropy with 
specific heat, volume expansion and compressibility. 
All these plots manifest an infinite divergence at the 
critical point. Furthermore, the smeared nature of the 
phase transition gets highlighted. 

Due to this smeared transition a direct application 
of the Gibbsian approach (Ehrenfest's scheme) to ex- 
hibit and classify phase transitions fails. We thus de- 
velop a generalized Ehrenfest scheme to account for 
this smearing. We show that such a generalization is 
compatible with our graphical analysis. 

The discovery that black holes are indeed thermo- 
dynamical systems is built on a mathematical analogy 
between laws of black hole mechanics and laws of ther- 
modynamics [2l. For chargeless, rotating black holes 
the first law of black hole mechanics relates the in- 
finitesimal change in black hole mass (Af) with the 
infinitesimal changes in its horizon area {A) and an- 
gular momentum (J), given by 



5M = —5A + 0(5 J, 
Stt 



(1) 



where k, and Q, are surface gravity and angular velocity 
at the event horizon, respectively. This is very much 



analogous to the first law of thermodynamics, 
dE = TdS - PdV. 



(2) 



The fact that the internal energy {E) of black holes 
being represented by its mass (M), together with 
Hawking's discovery that black holes have tempera- 
ture T = set the analogy between ([T]) and ([2]), 
associating an entropy to black holes S = ^. Several 
alternative approaches [3] have been developed which 
reproduce identical results for black hole temperature 
and entropy. 

We observe that while the first terms in the right 
hand sides of ([T]) and ([2 ) set a mathematical analogy 
between black noles and ideal fiuids, it is the second 
terms which actually establish the physical motiva- 
tion. It leads to the correspondence P < — > —VI and 
V i — > J. One can realise the physical importance 
of this analogy in the following manner. In ideal fiu- 
ids, both for isothermal and isentropic processes if one 
increases the pressure its volume shrinks so that the 
density increases. Thus the P — V plot for fluids has 
a negative slope. On the other hand the relation be- 
tween the angular velocity and angular momentum of 
Kerr-AdS black hole is given by [H [5] 



(7r+5')(7r+5-502) 



27r3/2(7r + S')3 



(3) 



We plot the fl ~ J graph in figure- 1 by taking vari- 
ous constant values for entropy. All these plots have 




FIG. 1: Angular velocity (Q) vs. angular momentum (J) 
plots. For blue line (lower )S = 1; for purple line (middle) 
S — 3; and for grey line (upper) 5* = 5. 

positive slopes which give a physical credibility to the 
correspondance, 
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(4) 
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FIG. 2: Entropy (5') and Gibb's free energy (G) plot with 
respect to temperature (T) for fixed $1 = 0.5 



Now we are in a position to use the known tools of 
thermodynamics to black holes. In order to follow a 
Gibbsian approach to exhibit and classify black hole 
phase transitions it is customary to define the Gibb's 
free energy. Gibb's free energy and temperature of 
this black hole was calculated earlier in [4ji5j and these 
are given by 



G{s, n) = M -TS - nj 

S{TT + Sf 



s-sn^ 



2TrS(n^ - 2) ~ 3S^{n^ - 1) 
47riS'(7r + S')2 



(5) 



(6) 



The relevant S — T and G — T plots are shown in fig. 2. 
The temperature-entropy plot is continuous. Also, 
there exists a minimum temperature (Tc — 0.2574) 
after which the slope of this curve tends to change its 
sign very slowly. These features are qualitatively sim- 
ilar to fluids (or magnetic systems) |S] revealing that 
there is no first order phase transition involving la- 
tent heat. Also, the same analogy indicates the onset 
of a continuous (higher order) phase transition with a 
critical temperature Tc- Contrary to what happens in 
usual fluids, however, the transition appears smeared 
instead of being sharp. Similar conclusions are also 
drawn from the G — T plot. 

The situation becomes more transparent when, 
once again following the analogous treatment for flu- 
ids, we plot the specific heat at constant angular ve- 
locity (Cn) (analogue of Gp for fluids) with entropy. 
The expression of Cn is given by |5j, 

Cn = 2S{tt + S){tt + S- Sif) x 



2TTS{n'^ - 2) - 352(172 - 1) 



(tt + SfiiS - tt) - 652(71 + Syn^ + 53(47r + 35)1^4 



Second, it is seen that the curvature of the two arms 
in fig. 2 (G — T) change sign. Since the curvature of 
G is just the specific heat, we observe the occurence 
of the phase transition curves in the positive and neg- 
ative quadrant of the Cn — S plot seperated by the 
critical point. 

Clearly a Gibbsian approach, which has its root 
in mean field theory, breaks down in order to analyse 
this infinitely diverging, smeared phase transition. In 
the following we make a modification to known Gibb- 
sian approach (Ehrenfest's scheme [7]), incorporating 
a smearing effect, to analyse this new phase transition. 
We shall build our machinary considering a fiuid sys- 
tem which undergoes a smeared phase transition and 
finally using Q we shall recast relevant quantities for 
black holes. 

Now to generalize Ehrenfest's scheme we start from 
the fact that in a smeared phase transition the values 
of entropy or volume of two phases at the phase tran- 
sition point are not exactly same as it should be in a 
second order phase transition, instead there are small 
smearing terms (s,v) in the following manner 



The Cn — S plot (figure 3) is a concrete evidence of on- 
set of a phase transition. Now if we want to compare 
this plot with usual higher order phase transition we 
find one important difference. This is the appearance 
of a smeared region on either side of the critical point. 

Let us just briefly discuss how the effects mani- 
fested in the curves of fig. 2 are concretised in fig. 3, 
which is quite akin to what happens for fluids, albeit 
with some crucial difference. First, the smeared na- 
ture of the phase transition is prominently displayed. 



S2 — Si 



(8a) 
(8b) 



where the subscripts 1 and 2 denote the values in the 
two phases. If we characterise the phase transition by 
the temperature (T) and pressure (P) for which ( [8a| ) 
and (8b I hold, at a different phase transition pomt 

dT, P' 



(characterised by T' = T 
following equations will be true 



S2 
V2 



dS2 
dV2 



Si 

Vi- 



-dSi 
dVi - 



- s 
v' 



P + dP) the 

(9a) 
(9b) 



From (8a I and (9a) we find 

dSi = dS2 + (s 



(10) 



Taking 5 as a function of T and P we write the in- 
finitesimal change in entropy as dS 



-.dT ) 



.dT 



'^^^^dP. Using the Maxwell relation (§|^)p 



•.dP )t 



this equation takes the form, 



dS 
where, 

X 



Cf 



dT - xdP, 



df 



Cp^T 



dS_ 
df 



(11) 



(12) 



are, respectively, the volume expansion index and the 
specific heat at constant pressure. Since (111 is inde- 
pendently true for two phases we write 



(7) 



dSi = ^dT-xidP; 



dS2^^dT^ 



X2dP.{13) 



Now substituting this in (10 1 we finally get the gen 



eralized first Ehrenfest's equation 

'dP\ ^Cp,-Gp, 
dT, 



T(x2-Xi) (X2-Xi)(r'-T)' 



(14) 



In the absence of any smearing in entropy, the second 
term on the right hand side is zero and one recovers 
the first Ehrenfest's equation in its usual form. 
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Next considering ( 8b ) and ( 9b I we find 



dVi ^dV2 + {v-v'). 
Taking F as a function of T and P, we find, 



dV 



dV 
df 



dT 



dV 
dP 



dP 



= xdT + ^dP 
where ^ is the compressibihty 



dV 
dP 



(15) 

(16) 
(17) 

(18) 



and X has been defined in (12). Following the previ- 
ous steps we find the generalized second Ehrenfest's 
equation 



dP 
df 



X2 - Xi 



{v - v') 



6-6 (6-6)(r'-T)- 



(19) 



We can e xpe ct that in a smeared phase transition, 
(14) and (19) will be satisfied. In the following part 
of our paper we shall investigate this in the context of 
black holephase transition. 

Using (Hi), ( [l4| ) and ( 19 ) we now propose the follow- 
ing generahzedEhrenfesFs equation for the Kerr-AdS 
black hole, 



df 
df 



_ - Cqi s - s' 

^Tix2-Xi) iX2-Xi)(T'-T) 

(20) 

^ X2-X1 a - f) ^ 

6-6 (6-6)(T'-r)- ^ ' 

In order to check the validity of ( [20| and ( [2l] ) we 
first recall the following equations from [5j 

X = 

652(7r + Sfn - 2S'3(7r + S')(27r + 3S)n^ 
(tt + S')3(35 - vr) - 652(7r + S)^n^ + S^4tt + 3S')f74 

(22a) 



(22b) 



an 

df 

dn 

df 



4tt2 {tt + s - sn^ 



2^2 / s{^+sr 



Sn iS{TT + 5')(27r + 35)172 _ 3(7^ + S)^) 

(22c) 





FIG. 3: Semi-classical Specific heat (Cn) and volume ex- 
pansion index (x) vs entropy (S) plot. The solid line 
stands for Q, = 0.5 and the dashed line corresponds to 
its infintely small variation (Q' — Q + 5% x Q). 



FIG. 4: Moment of inertia index (5) and angular momen- 
tum (J) vs entropy (5*) plot. The solid line stands for 
$1 = 0.5 and the dashed line corresponds to its infintely 
small variation {Q' — Q + 5% x Q). 



(Sc = 1.208). For these values of 12 and S the critical 
temperature is found to be Tc — 0.2574 from 

Now t he le ft hand sides of ( |20| , {2l\ evaluated 
by using (22c I, (22d| at the critical 
Tot n = 0.5) 



using 
1.208, 



point {Sc 

U'2570or n = 0.5) are found to be 
12.8027 and 12.7925. To calculate the right hand side 
we first draw the Cn — S , \ — S . $— S and J— S plots 
(figures 3,4) using ([22a| , ( |22b[ ) and respec- 
tively. Le t us fi rst consider the smearing independent 
terms in ( 20 ) , ( [2T| ) . The respective values of relevant 
quantities are borrowed from their plots. We choose 
such points in two phases (say 5*1 for phase 1 and 
iS'2 for phase 2) by avoiding the smearing region as 
in this part values widely change for two neighbour- 
ing points. Our prescription of choosing a point is 
that it should be very close to the critical point and 
belong to a region where the slope of the curve is van- 
ishingly small. For — 0.5 we take Si — 0.653 and 
5*2 — 1.763. Thus we hav e a smeared region of width 
s — S2 — Si — 1.11 (see (8a)). For this set of values 



in + s-sn^ 

5(7r + 5)3 



H 



(22d) 



the changes in Ca,x,^ between the two phases are 
found to be ACq = 23.80 , Ax = 7.72, = 0.90. 
The values thus found for the first terms of the right 
hand side of pO]) and ^ are 11.9817 and 8.578 re- 
spectively. Wenow observe that, since the l.h.s. of 



(20) equals 12.8027, the r.h.s. almost matches. This 



where F and H are defined as P]. The physically al- 
lowed values for are < < 1, else the Hawking 
temperature becomes negative 0. We choose a par- 
ticular value around the central range (fi = 0.5) to 
examine how this generalized scheme works. For this 
value of n the plots of Cn, X ^'^d ^ (see figures 3 and 
4) show a smeared divergence at the critical entropy 



does not hold for (21). This is reminiscent of recent 
studies fS] which show that if one neglects the smear- 
ing effect for black hole phase transitions then only 
the first Ehrenfest's equation is satisfied and the sec- 
ond one is violated. Effectively, therefore the smearing 
terms have a nontrivial contribution only for (21 ). We 
now discuss this issue. 

In order to obtain the second (smearing depen- 



n 


1st generalized Ehrenfest's rein. 


2nd generalized Ehrenfest's rein. 




l.h.s=-(|^)^ 






r.li.s=Xi + Yi 


l.h.s=-(|y)g 




Y2 


r.h.s=X2 + Y2 


0.1 


72.1469 


64.6853 


0.0 


64.6853 


66.8405 


6.232 


39.50 


45.732 


0.2 


35.6305 


31.6156 


Do 


31.6156 


35.5666 


9.470 


20.56 


30.030 


0.3 


23.2151 


20.5940 


Do 


20.5940 


23.2601 


10.329 


12.36 


22.689 


0.4 


16.8036 


15.0467 


Do 


15.0467 


16.7774 


9.833 


5.13 


14.963 


0.5 


12.8027 


11.9817 


Do 


11.9817 


12.7985 


8.578 


2.924 


11.502 


0.6 


9.9750 


8.9837 


Do 


8.9837 


9.9852 


7.484 


1.26 


8.744 


0.7 


7.7547 


7.1375 


Do 


7.1375 


7.7544 


5.944 


0.596 


6.540 


0.8 


5.8653 


5.4616 


Do 


5.4616 


5.8674 


5.025 


0.262 


5.287 


0.9 


4.0034 


3.3300 


Do 


3.3300 


4.0030 


2.970 


0.188 


3.158 



TABLE I: Summary of results for validity of generalized Ehrenfest's relations for Kerr-AdS black hole for various angular 
velocities. 



dent) terms we first make a small change (5%) in an- 
gular velocity and see tlie behaviour of all the plots. 
Let us first consider ( 20 ). It is clear from all plots that 
the smearing width suffers no significant change, i.e. 
s — s' sa 0. We may neglect this small change as it is 
further weakend by the denominator of {s — s') in (20 1. 



ernment of India, for financial help. 



As a result the right hand side of the first equation al- 
most remains unchanged. Regarding (21), however, 
this small change in f2 changes the smearing in J. It 
is obtained by finding j{= J{S2 = 1.763,0 = 0.5) — 
J{Si = 0.653, f7 = 0.5)) and /(= J{S2 = 1.763, = 
0.525) - J{Si = 0.653, n = 0.525)) from ^ and then 
taking their difference. This value of j — j in ( plj ) is 
comparable with its denominator. Thus it is no longer 
negligible. The contribution from this extra term is 
found to be tc — /T/t., ^p^ = (o.09i o.ose 
This value when added with 



gives the right hand side of (21 ) as 11.502 
reasonable agreement with 
In order to get a full picture we consider va riou s al 



— 2 Q24 

0.0019x0.9 ^.c^-i. 

the other term (8.578) 
This is in 
tiel.h.s. value (12.7925). 



lowed values of and check the validity of (20) and 
(21). The results are summarised in a tabular form 
(see Table-I ^lOj). It clearly shows the validity of these 
new equations for black hole phase transition. 

We have discussed a new approach to study phase 
transitions in black holes that is based on their anal- 
ogy with fluids. Several features of this analogy played 
a crucial role. Unlike the liquid to vapour transition, 
here it is not first order. Rather, according to our 
analysis this is a smeared continuous (higher order) 
phase transition. Following a generalized Gibbsian 
approach we have derived the governing equations for 
such transitions. These equations were verified for 
phase transition in Kerr-AdS black hole. Our calcu- 
lations are robust. We checked this by taking other 
values of Si and S2 close to those taken here. We 
feel that the approach discussed here can be pushed 
further, by exploiting the analogy with fluids, to 
gain a deeper insight into black hole thermodynamics. 
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Jgj p — (T-3S)(7r + S)3-2S(7r + S)^(47r + 3S)n^+S^(2ir + 3S)^n* 

27r 2 [(7r-3S)(7r+S)* + 6S2(7r + S)3n2-S3(7r + S){47r+3S)n<l] 
_ 47r3/2s(,r+S)^Q((3(7r + S)2-S(27r + 3S)n2)) 

~ (7r + S)3(3S-7r) + 2S(7r + S)2{47r + 3S)Q2_s2(2^+35)2Q4 

[10] The new variables in table-II are defined as Xi = 

'-^'^2 — s-s' y _ X2-X1 yi _ 
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